\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{Sech}^{2}$$\end{document}Sech2-type solitary waves and the stability analysis for the KdV–mKdV equation

In this article, we investigated the solitary wave solutions of the KdV–mKdV equation using Hirota’s bilinear method. Closed-form analytical single and multiple solitary wave solutions were obtained. Through qualitative methods and the analysis of solitary waveforms, we discovered that in addition to sech-type solitary waves, the system also contains \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{Sech}^{2}$$\end{document}Sech2-type solitary waves. By employing the trial functions method, we obtained a single \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{Sech}^{2}$$\end{document}Sech2-type solitary wave and verified its existence and stability using the split-Step Fourier Transform method. Furthermore, we use the collision of two \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{Sech}^{2}$$\end{document}Sech2-type single solitary waves to excite a stable \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{Sech}^{2}$$\end{document}Sech2-type double solitary wave. Similarly, we excite a stable triple solitary wave with three \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{Sech}^{2}$$\end{document}Sech2-type single solitary waves. This method can also be used to excite stable multiple solitary waves. It is shown that these solitary wave solutions enrich the dynamic behavior of the KdV–mKdV equation and provide methods for solving \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{Sech}^{2}$$\end{document}Sech2-type solitary waves, which hold significant theoretical value.

In this article, we investigated the solitary wave solutions of the KdV-mKdV equation using Hirota's bilinear method.Closed-form analytical single and multiple solitary wave solutions were obtained.Through qualitative methods and the analysis of solitary waveforms, we discovered that in addition to sech-type solitary waves, the system also contains Sech 2 -type solitary waves.By employing the trial functions method, we obtained a single Sech 2 -type solitary wave and verified its existence and stability using the split-Step Fourier Transform method.Furthermore, we use the collision of two Sech 2 -type single solitary waves to excite a stable Sech 2 -type double solitary wave.Similarly, we excite a stable triple solitary wave with three Sech 2 -type single solitary waves.This method can also be used to excite stable multiple solitary waves.It is shown that these solitary wave solutions enrich the dynamic behavior of the KdV-mKdV equation and provide methods for solving Sech 2 -type solitary waves, which hold significant theoretical value.
The study of solutions for nonlinear partial differential equations (NLPDEs) is essential to the theory of physical mathematics and theoretical physics 1 .These solutions help us to better understand the mechanisms of complicated nonlinear physical phenomena and the dynamic processes modeled by these nonlinear evolution equations 2 .The KdV equation (Korteweg-de Vries equation) and the mKdV equation (modified Korteweg-de Vries equation) are two important types of nonlinear partial differential equations, widely applied in many physical systems.In some cases, these two equations can be combined to better describe complex physical phenomena.The combined KdV-mKdV equation is an important tool for describing nonlinear wave phenomena, with widespread applications in many fields, such as solid-state physics, plasma physics, fluid physics and quantum field theory [3][4][5] .By studying these equations, we can gain a deeper understanding and predict the nonlinear behavior in complex systems, thereby advancing related science and technology.The combined equations not only enrich the theoretical foundation of nonlinear science but also provide powerful tools for solving practical problems.These combined equations play a significant role in the study of nonlinear waves, solitons, and other complex systems.
The existence of soliton-type solutions for nonlinear models is highly significant due to their potential applications in various physics areas, including nonlinear optics, plasmas, fluid mechanics, condensed matter, and many more.In fact, many kinds of soliton solutions have been obtained by using for example, the inverse scattering method 6 , the Sub-ODE method 7 , the tanh method 8 , the Darboux transformation 9 , the Lie group method 10 , the algebrogeometric method 11 , the homogeneous balance method 12 and the Hirota's method 13 , The extended rational expansion method 14 , The MSE method 15 and so on.
The combined form of the Korteweg-de Vries equation and the Modified Korteweg-de Vries equation form the following so-called combined KdV-mKdV equation 3,1 : where α , β and γ are constants.Specifically, α controls the strength of the first-order nonlinear effects in the equation, β adjusts the strength of the second-order nonlinear effects, typically describing stronger nonlinear effects, γ controls the strength of the dispersion effects of waves.
The KdV-mKdV equation is completely integrable, giving rise to multiple-soliton solutions.The combined KdV-mKdV equation and its generalizations have attracted increasing interest, as evidenced by numerous studies [16][17][18][19][20][21][22] .The authors obtained tanh-sech type solitary solutions and triangular functions solutions in (1) www.nature.com/scientificreports/KdV-mKdV equation in Ref. 23 .In Ref. 24 , the tanh, coth and tan type solitary wave solutions were studied.In Ref. 5 , the periodic progressive wave, and one and two soliton solutions were studied by various analytical methods.In Ref. 25 , the bell-shaped soliton and a new soliton solution were obtained by a kinds of Jacobi doubly periodic solutions.Moreover, new periodic wave solutions were constructed by a method with the aid of a Sub-ODE.The rational solutions with free multi-parameters and N-soliton solutions were obtained by the bilinear approach in Ref. 26 .Although nonlinear solitary waves in the KdV-mKdV equations have been studied extensively both analytically and numerically, there are still many open problems.For instance, we have well-developed methods for solving and analyzing soliton waves of the sech type.However, to our knowledge, there is no research on Sech 2 -type solitary waves, including the existence and stability of Sech 2 -type solitary waves in the KdV-mKdV equations.
In this work, we consider the generalized combined KdV-mKdV equation (1).By employing the method of trial functions, we derived several approximate analytic solutions of Sech 2 -type solitary waves.Through the application of perturbation methods 27,28 , we investigated the stability of both single and multiple Sech 2 -type solitary waves, and successfully obtained stable solitary wave solutions for both single and multiple Sech 2 -types.

The bilinear form
The Hirota's bilinear method is an effective approach for obtaining analytical solutions to nonlinear evolution equations.Therefore, we used Hirota's bilinear method to solve Eq. (1).
In order to find the soliton-like solutions of Eq. ( 1) we set where a 0 is a constant to be determined later, and Substituting the Eq. ( 2) into Eq.( 1) yields Integrating both sides of Eq. ( 3) by x once and setting the constant of integration equal to zero yields: Expanding the Eq. ( 4) leads to By a simple transformation, Eq. ( 5) can be transformed into the following equation: By substitution Eq. ( 6) into Eq.( 1) and using the properties of the bilinear operator, we obtain In order to obtain the bilinear form of Eq. ( 7), setting β = 6 , then we have Further, we obtain the bilinear form of Eq. ( 8): It should be noted that the bilinear form Eq. ( 9) is obtained when the coefficient β = 6.
According to the Hirota's bilinear method , we may choose the f and f * in the form (5) www.nature.com/scientificreports/Substituting Eq. ( 10) into Eq.( 9) and collecting coefficients of polynomials of with the aid of Mathematica, then equating the coefficients of various powers of ε to zero, we get a set of algebraic equations: Suppose that where ξ 1 = ω 1 t + k 1' x + c 1 , θ 1 are constants to be determined later Substituting Eq. ( 14) into Eq.( 11), we have

Analytical solution
Based on the bilinear form of the KdV-mKdV equation, we can obtain analytical solutions for single and double solitary waves.

Single solitary wave solutions
For single soliton solutions we set = 1 + e ξ 1 +iθ 1 and ε = 1 , then we have Substituting Eq. ( 17) into Eq.( 2), the single soliton solutions are given by which can be rewritten as where

Double solitary wave solutions
Suppose that where ξ i = ω i t + k i' x + c i , θ i , i ∈ {1, 2} are constants to be determined later.Substituting Eq. ( 21) into Eq.( 12), we have For double soliton solutions we set (1) .

Sech 2 -type solitary wave solutions
Using 0 ≤ |cos θ | ≤ 1 from Eq. ( 19), we have and it can be transformed into From the above qualitative analysis, it is easy to conclude that the waveform of solution (19) lies between that of sech-type and Sech 2 -type solitary waves.
When sin θ 1 = 1 , solution (19) reduces to the standard sech-type solitary wave solutions In the limit case when θ 1 → 0 ( sin θ 1 ∼ θ 1 , cos θ 1 ≈ 1 ) the solution (1) become the Sech 2 -type solitary wave solutions It is worth noting that due to θ 1 possibly approaching zero, the solution may have significant errors.Solution (28) represents a specific analytical solution.To obtain a more general Sech 2 -type solitary wave solution, we need to adjust the parameters of solution (19) and solve again.
Suppose that which can be rewritten as where , a 0 , a and M are constants to be determined.Substituting Eq. ( 30) into Eq.( 1), we obtain Further simplification of Eq. (31) yields: Vol Collecting all terms with the powers in s ech −n ξ(n = 0, 1, 2, 3) and setting each of the obtained coeffificients for s ech −n ξ to zero yields the following set of algebraic equations with respect to a 0 , k 1 , ω 1 , α, β, θ 1 and M : Solving Eqs. ( 33)-(36), we obtain the following solutions: Case (i): The travelling wave solution of KdV-MKdV equation are given by (a) When θ 1 = π/2 , then where (37) (39) where where

The linear stability of solitary wave solutions of Sech 2 -type
The Split-Step Fourier Transform (SSFT) method is a pseudospectral numerical technique utilized for solving nonlinear partial differential equations.This method involves computing the solution in incremental steps.It handles dispersive and nonlinear effects independently through a sequential approach: first addressing the nonlinear effects and then the dispersive effects.
In this section, we will use the split-step Fourier transform (SSFT) method 27 to study the linear stability of single and multiple Sech 2 -type solitary waves in equation (1).To this end, we use the waveform of a Sech 2 -type solitary wave (refer to solution (40) or (44)) at time t=0 (selecting waveforms at other times is also acceptable and does not affect the analysis results) with a perturbation of a random uniformly distributed noise field of amplitude 10 −4 as the initial condition.By observing the changes in perturbations over time, we can determine the stability of the solitary wave.If the perturbations are suppressed, the solitary wave is stable; if the perturbations increase exponentially, the solitary wave is unstable [27][28][29][30] .Solution ( 28) is obtained under θ 1 → 0 , but a key issue is that the wave number k 1 must be very large to ensure the validity of the solution.Therefore, we solved Equation (1) again and obtained Solution (40) or (44), which ensures a reasonable value for k 1 .
Based on the above reasons, it is reasonable to use Sech 2 -type solitary waves as the input signal without con- sidering the dispersion relationship to study the stability of solitary waves.
Figure 1a shows a stable train of solitary waves excited by a Sech 2 -type solitary wave as the initial condition, with nearly thirty solitary waves from left to right, gradually decreasing in amplitude.From the perspective of energy conservation, the waveform of the solitary wave with the maximum amplitude (on the far left) no longer maintains a Sech 2 shape.Therefore, although we used a Sech 2 -type solitary wave as the initial condition, stable Sech 2 -type solitary waves were not excited.
Relative to Fig. 1a, the solitary waves in the stable train shown in Fig. 1b not only exhibit reduced amplitude variations but also decrease in number.The amplitude of the excited solitary wave (on the far left, with the maximum amplitude) in Fig. 1b more closely resembles a Sech 2 waveform.
The stable solitary wave excited in Fig. 1c does not show any changes compared to the initial waveform, clearly indicating the excitation of a Sech 2 -type solitary wave.
Relative to Fig. 1a,b, the stable solitary wave excited in Fig. 1d (on the far left) has a larger amplitude, and the train of solitary waves in Fig. 1d contains fewer solitary waves with smaller amplitudes.It is evident that the solitary wave excited in Fig. 1d does not exhibit a standard Sech 2 shape.However, it is a stable solitary wave.
In Fig. 1e, we can see the initial solitary wave splitting into two solitary waves.Before t = 3.5 , both of these solitary waves maintain stable propagation.However, at t = 3.5 , the solitary wave on the left suddenly increases in amplitude and changes direction.It is evident from the figure that stable solitary waves were not excited.
In Eq. ( 1), the amplitude of the solitary wave is a factor that affects the stability as shown in Fig. 1.It can be observed that for the same solitary wave, as the wave amplitude increases gradually, the solitary wave becomes stable initially.However, with further continuous increase in wave amplitude, the solitary wave begins to become unstable.
We use the conclusion from Ref. 29 , "the double-soliton is inspired by the elastic collision of two singlesolitons", to construct a Sech 2 -type double solitary wave.
In Figure 2a shows stable Sech 2 -type double-solitary waves excited by the interaction of the initial single solitary waves u D 1 and u D 2 .It is evident that the amplitudes of u D 1 and u D 2 do not change, indicating that the excited double solitary are of the Sech 2 type.Figure 2b utilizes the same method as Fig. 2a, using the interaction of single solitary waves u D 1 , u D 2 , and u D 3 to excite triple solitary waves.It is easy to see from the figure that stable triple solitary waves of the Sech 2 -type are excited.Figure 3 shows the stable quadruple solitary and quintuple solitary waves excited using the same method as in Fig. 2. Of course, if the initial conditions are chosen improperly, unstable multiple solitary waves may also be excited.

Conclusion
This paper investigates the solitary wave solutions of the nonlinear KdV-MKdV equation using the Hirota's bilinear method, focusing on both single solitary wave solutions and multiple solitary wave solutions.We know that the nonlinear terms of the KdV and mKdV equations often coexist in certain physical systems, such as fluid dynamics and quantum field theory, giving rise to the KdV-mKdV equation.Interestingly, in the general KdV equation, exact solitary wave solutions of the Sech 2 -type exist, whereas in the mKdV equation, exact solitary wave solutions of the sech-type are present.The KdV-mKdV equation supports exact solitary wave solutions of the sech-type but does not support general exact solitary wave solutions of the Sech 2 -type.Therefore, investigating the existence and stability of solitary waves of the Sech 2 -type, both single and multiple, is of significant importance.
In the limit of analytic solutions, the Sech 2 -type solitary waves can be obtained.However, this limita- tion requires a very large wave number.To obtain general Sech 2 -type solitary waves, this paper re-solves the KdV-mKdV equation using the method of trial functions, resulting in general Sech 2 -type solitary waves.Addi- tionally, the propagation stability of single and multiple Sech 2 -type solitary waves is analyzed and simulated using the SSFT, leading to stable single and multiple solitary waves.

Figure 2 .
Figure 2. Stability analysis of Sech 2 -type double-soliton or triple-soliton inspired by a perturbed excitation with α = −6, β = 6, γ = 4.887 : (a) initial condition selected from two single-soliton solutions u D 1 and u D 2 ; (b) initial condition selected from three single-soliton solutions u D 1 , u D 2 and u D 3 .